Theory of the high-frequency chiral optical response in a p x + ip y superconductor 
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The optical Hall conductivity and the polar Kerr angle are calculated as functions of temperature 
for a two-dimensional chiral p x + ip y superconductor, where the time-reversal symmetry is sponta- 
neously broken. The theoretical estimate for the polar Kerr angle agrees by the order of magnitude 
with the recent experimental measurement in S^RuCU by Xia et al. [PRL 97, 167002 (2006)]. The 
theory predicts that the Kerr angle is proportional to the square of the superconducting energy gap 
and is inversely proportional to the cube of frequency, which can be verified experimentally. 

PACS numbers: 74.70.Pq, 78.20.Ls, 74.25.Nf, 73.43.Cd 



Xia et al. [l[ recently reported experimental observa- 
tion of the polar Kerr effect in the superconducting state 
of Sr2RuC>4. In the absence of an external magnetic 
field, reflected light shows rotation of polarization, which 
is a clear signature of the spontaneous time-reversal- 
symmetry breaking in the superconducting state Q ■ Pre- 
vious muon spin relaxation measurements Q suggested 
the time-reversal symmetry is broken in Sr2Ru04, but 
the polar Kerr experiment [lj gives a much more con- 
vincing evidence for this remarkable effect. 

Sr2RuC>4 consists of weakly coupled two-dimensional 
(2D) metallic layers. It was proposed theoretically that 
the superconducting pairing in this material is spin- 
triplet Q and has the chiral p x + ip y symmetry [Bj . Such 
an order parameter breaks the time-reversal symmetry 
and is analogous to the 2D superfluid 3 He-A [6j. There 
is substantial experimental evidence in favor of the spin 
triplet and odd orbital symmetry of the superconducting 
pairing in Sr2Ru04 Q, which includes measurements of 
the spin susceptibility Q and the Josephson effect Q 
(see, however, an alternative interpretation [lfj). On the 
other hand, the chiral character was not so well estab- 
lished experimentally (see Ref. [ljj for interpretation of 
tunneling measurements). 

Although the experimental demonstration |l| of the 
spontaneous polar Kerr effect is very convincing, a the- 
ory of this effect for chiral superconductors is not well 
developed. Theories 12, Tj^ concluded that there is no 
chiral term in the single-particle response of a p x + ip y 
superconductor, although Fig. 1 of Ref. [l2| shows a non- 
zero Kerr effect for a different state, and Ref. [k| found 
some chiral response from collective excitation. On the 
other hand, Ref. 0] obtained the intrinsic quantum Hall 
effect in the single-particle response of a p x + ip y super- 
conductor, which was then studied in much detail in Ref. 
11]. Following Refs. HE EH, the polar Kerr angle 9 K 
can be expressed in terms of the imaginary part of the 
ac Hall conductivity a" v (Q) at a frequency f2 
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in terms of the 2D Hall conductivity a xy per one layer, 
which is related to the 3D one via the interlayer distance 
d. The natural dimensional scale for the 2D a xy is e 2 /h. 

In this paper, we calculate the ac Hall conductivity 
a xy (Q) at a finite frequency as a function of tempera- 
ture T for & p x + i p v s uperconductor. We generalize the 
method of Refs. @,13] and obtain the Chern-Simons-like 
term in the effective action at finite T and f2. In the inter- 
mediate calculations, we set the Planck constant to unity 
ft — ► 1, but restore it in the final results. The Lagrangian 
of electrons L = idt — H (where H is the Hamiltonian) 
for the 2D p x + ip y superconductor has the form [6j 



L 



id t + V 2 /2m + «(V- 
i(V •**+**• V)/2 id t 



* + * • V)/2 
- V 2 /2to - fi 



(2) 

Here d t and V = (d x ,d y ) represent time and space 
derivatives, m and fi are the mass and the chemical po- 
tential of the electrons. We assume the parabolic disper- 
sion law e(p) = p 2 /2m — fjb, where p = (p x ,p y ) is the 
electron momentum. The superconducting order param- 
eter is * = A x x + iA y y, where x and y are the unit 
vectors in the x and y directions. Because of the square 
symmetry in Sr2Ru04, we have A x = A y , but it is con- 
venient to label the two components of differently for 
the clarity of calculations. In momentum representation, 
Eq. ^ can be written as 



L = iu — e(p)r 3 - p x A x T\ +p y A y T 2 , 



(3) 



n(n 2 - l)Cld" xy 
where n is the refraction coefficient. We write Eq. ([I]) 



where the Pauli matrices r act on the spinor 
[ip(p), ip + (—p)] consisting of the particle and hole op- 
erators. We do not write the spin indices of electrons 
explicitly. The two spin components give the same con- 
tributions to the Hall conductivity, so the final results 
should be multiplied by 2. However, by introducing elec- 
tron and hole operators, we artificially doubled the num- 
ber of components, so the final result should be divided 
by 2 0]. Thus, we can obtain the correct result by consid- 
ering just one spin component, as implied in Eq. (|3|). In 
Eq. ((3|), we use the Matsubara frequency iu>, because we 
will be doing calculations at a finite temperature. From 



2 



Eq. ([3]), the Green function of electrons G = L 1 is: 



G(p) 



iu) + e{p)r 3 + p x A x n - p y A y r 2 
lu 2 + E 2 (p) 



(4) 



where p = (uj,p) is the three-component frequency 
momentum vector, and E{p) 



le 2 (p)+p 2 x A 2 x +p 2 y A 2 y 

is the electron dispersion in the superconducting state. 

To calculate electromagnetic response of the sys- 
tem, we introduce the electromagnetic potentials A = 
(A , A x , A y ) by using the long derivatives — iV =F eA/c 
and —idt ± eAo in the diagonal terms of Eq. @, where 
e is the electron charge, and c is the speed of light Q. 
We also assume that the superconducting order param- 
eter has a space-time-dependent phase so that it can 
be written as *f? = e^V&o, where is uniform with the 
real A x and A y . We will see that the effective action 
depends only on gradients of p. To the first order in A 
and ip, we find the following addition to the Lagrangian: 



r = 



-eA - i 
V-^p + p 



eA-V 



V ■ * <p + V *o • V 



V 



lAn - 



.eA-V 



(5) 



The Fourier transform of Eq. ([5]) can be written as 

T = -eA Q r 3 + eA ■ p/mc - p>p x A x r 2 - pp y A y T X . (6) 

Here the variables A and ip are assumed to be functions 
of the Fourier variable q = (Q, q x , q y ). Thus, the vertex 
T(q,p) ([6]) is a function of two vector arguments. 

The effective action of the system to the second order 
in r is 



Substituting Eqs. ^ and ^ into Eq. ([7]), we write S in 
the form 



Ci 

c 2 



(8) 



q.p 



where the denominator is 

C 2 = [{uj+n/2) 2 +E 2 {p+q/2)][{u-n/2) 2 +E 2 {p- q r2)\ 
and the numerator is 



(9) 



C\ = - Tr [~eA (q^T 3 + ep x A x (^/mc + ep y Ay(q^/mc-p(q)p x A x T 2 -p{q)p y A y T 1 } 

x [i(to + Q/2) + e(p + q/2)r 3 + (p x + q x /2)A x r 1 - ( Py + q y /2)A y T 2 ] 

x [-eA (-q)T 3 + ep x A x (-q)/mc + ep y A y (-q)/mc- tp(-q)p x A x T 2 - p(-q)p y A y Ti] 

x [i(uj-n/2) + e{p-q/2)T 3 + {p x -q x /2)A x T 1 -{p y -q y /2)A y T 2 ] (10) 

I 



The calculation of the effective action ([7]) is conceptually 
similar to the calculation of electromagnetic response in 
the BCS theory of superconductivity [12|,|l6|]. However, 
we focus only on obtaining the Chern-Simons-like term 
responsible for a xy @, Picking the Aq term from 

the first factor in Eq. (|10p and the A x or A y term from 
the third factor, we obtain a non-zero contribution after 
taking trace over the r matrices. The same procedure 
works for the A x or A y term from the first factor and the 
Aq term from the third factor. Combining these terms 
and changing the variable of integration q — > — q in the 
latter term, we obtain one contribution to C\ 

C[ a) = A (q) [-q yAx {-q)pl + q x A y (-q)p 2 y ] 2l ^v e \ 

M mc 

(11) 

In deriving we omitted the terms proportional to 

the product p x p y , which would vanish after integration 
over p x and p y . The integration over momentum p in 
Eq. ([5]) is concentrated near the Fermi surface, so we can 



replace p 2 — > p 2 F /2 and p 2 — > p 2 F /2 in Eq. pT|) . because 
P x + Py ~ Pf- Making the Fourier transform of Eq. (fTTj) 
to the coordinate space, we find 

C[ a] = A (d y A x - d x A y )^, (12) 

mc 

where A = A x pp = A y pp is the energy gap at the 
Fermi level. 

Picking the last two terms in the first factor in Eq. 
(fTU)l and the A x and A y terms in the third factor or vice 
versa, we obtain another contribution to C%: 

C[ b) = tnp(q)[-q y A x (-q)p 2 x + q x A y (-q)p 2 y }^^ . 

" mc 

(13) 

Replacing the Matsubara frequency by the real frequency 
ifl — > Q in Eq. (|13p and Fourier-transforming to the co- 
ordinate space, we find 

C[ b) = d t p(d y A x - d x A y )^l, (14) 
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Combining the contributions (p~2|) and (fT4| to C% and 
substituting into l[8]). we find a Chern-Simons-like term 
in the effective action Ki LLJI 



5*cs = Oxy I dtdxdy(A + d t ip/2e)(d y A x - d x A y )/c, 

(15) 



where 



A^ 2 



'xy 



y- 



p 



(16) 



is the effective Hall conductivity. Indeed, taking the vari- 
ational derivative of Eq. (I15|) , we find electric current 



SS C s 



■nj 



E - ±d t ( V<p - —A 

2e \ c 



x z 



, (17) 



where E = —'VAq — dtA/c is the electric field, and the 
last term in Eq. (|17[) is proportional to the time deriva- 
tive of the London supercurrent j s = (p s e/2?7j)[V</? — 
(2e/c)A\. Obtaining a self-consistent equation of motion 
for the superconducting phase tp is a complicated prob- 
lem jl7| . However, one may argue that the supercurrent 
contribution in Eq. Q17p is ineffective at high frequencies, 
so the last term can be omitted, and we obtain the stan- 
dard relation for the Hall conductivity j = a xy E x z. 
The Chern-Simons-like term (|15|) was derived in Rcf. [6] 
at T = and in Ref. [3] near T c via the Ginzburg- 
Landau expansion. Notice that is does not have the com- 
ponent A x dtA y , which is present in the standard Chern- 
Simons term e /ilyA A^F V \. Nevertheless, Eq. (fT5|) is gauge- 
invariant, because the gauge transformation of A in the 
first factor is compensated by transformation of the su- 
perconducting phase ip, and the last factor is manifestly 
gauge-invariant Q. 

Now we substitute Eq. © into Eq. (fTE)) and derive an 
explicit expression for the Hall conductivity. To obtain 
optical response at a finite temperature, we do analyti- 
cal continuation from the Matsubara to real frequencies, 
which is well known in the BCS theory [l6j |. We take the 
limit q — > while keeping finite frequency, as appropriate 
for optical response, and find 



a xy {VL) 



ds 



1 



1 - 2n(E/T) 
E 2 



(18) 



1 



fl + Z7 - 2E tt + i~/ + 2E 



where E = y/ 'e 2 + A 2 ,, n(E/T) is the Fermi distribution 
function, f2 is the real frequency, 7 is a relaxation rate, 
and we replaced the integration over dp x dp y / (2ir) 2 m by 
the integration over de/2-K. First we take the dc limit 
f2 — > at zero temperature T — > in Eq. (TT8|) and find 



dc 
xy 



e 

An 



e 

2h' 



(19) 



where we restored the dimensional factor ft — h/2ir in 
the denominator. Eq. (fT9| demonstrates the half- integer 
quantum Hall conductivity in agreement with Ref. [g]. 
As discussed in Ref. [IH , it is difficult to measure the dc 
Hall conductivity experimentally because of screening by 
supercurrents. 

At a finite temperature T, the dc Hall conductivity 
Vxy( T ) = (e 2 /2h)f d {T) is reduced by the factor 



A 2 roo 

/-CO = ^ / de 
L J -00 



1 - 2n(E/T) 
W 2 ' 



(20) 



The factor fd{T) interpolates between 1 at T = and 
at T = T c and behaves as f d (T) a A a ^T c - T 
near T c . The same factor ([20)) describes temperature de- 
pendence of the quantum Hall effect in the magnetic- 
field-induced spin-density-wave (FISDW) state of the 
quasi-one-dimensional organic conductors (TMTSF^X 
flil Hit ]. Eq. (|20|) represents the dynamic limit of the 
dc electromagnetic response 

[IE If the limit n 
is taken in Eq. (TTI)]) first and then q — > 0, that would 
generate the static limit f s (T) for the reduction func- 
tion, which has the same temperature dependence as the 
London superfluid density p s (T), particularly f s (T) oc 
A 2 oc (T c — T) near T c (see discussion in Sec. VI of Ref. 
(HI). Ref. E 

by doing the Ginzburg-Landau expansion. 

Now we calculate the imaginary part of the Hall con- 
ductivity a xy (Q) at a high frequency ft ^> A . One con- 
tribution originates from the pole at £1 — 2E in Eq. (|18|) . 
This term represents creation of an electron pair above 
the energy gap Ao or a hole pair below the gap by absorp- 
tion of a photon with the frequency Q. By integrating 
over e « ±E in the vicinity of the resonance, we find 
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obtained the static limit for cr~; near T c 



'xy 



(fi) = 



2n 2 



e 
2ft 



Aq 
Tift 



(21) 



Calculating this term, we set E 2 = (ft/2) 2 in the denom- 
inator of the first factor in Eq. (fT5|) and put n(E/T) = 0, 
because ft 3> T. We observe that Eq. (f2"Tj) does not de- 
pend on the relaxation rate 7 and is reduced relative to 
Eq. |19|) by the factor (A /hil) 2 /2tt. The temperature 
dependence of a xy (Q) is given by Aq(T). 

There is another contribution to the integral (fT8|) orig- 
inating from the peak in the density of states at E « Ao. 
By changing the variable on integration from e to E, we 
rewrite Eq. (fl"8"|) as follows 



e 2 A 2 

Czy^) = 



dE 

A,, ^ 



1 - 2n{E/T) 



1 



a 2 (ft 



4E 2 
(22) 

Integral (f2"2"| over <il? is logarithmic between Ao and 
ft. For simplicity, we consider low temperatures, where 
ii{E /T) w 0, and find the following contribution a xy to 
the Hall conductivity 

2 A 2 / n 



<T xy (fl) 



A 2 



7T (fi + i7 ) 2 ^Aq 
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a xy 



(Q) 



4e 2 A 2 ft 7 
~~h (hCl) 3 



In 



(23) 



Eq. (|2"3"|) is reduced relative to Eq. (j2"Tj) by the factor 
and is enhanced by the factor ln(?iil/Ao). Using 
the numbers from Ref. [l|] and given below, we conclude 
that the reduction of Eq. (|23|) is much greater than the 
enhancement, so we focus only on Eq. (|2Tj) . 

Substituting Eq. (HTJ) into Eq. |T]), we find the Kerr 
angle 



2tt 



>K 



A 2 
^0 



A 



A 2 
^0 



n(n 2 ~l) d (Ml) 3 n(n 2 -l) d {HQ) 2 ' 



(24) 



where a = e 2 /tic — 1/137 is the fine structure constant, 
and A is the wavelength of light. Using the interlayer 
distance d = 1.3 nm [7j and the values n(n 2 — 1) = 3 
and A = 1550 nm from Ref. [l|, we find that the first 
two factors in Eq. give 2.9. Using the BCS for- 

mula A = 1.76 k B T c = 0.23 meV for T c = 1.5 K | 
and UVt — hc/X — 0.8 eV, we find that the last factor 
in Eq. (gj) is (A /hn) 2 = 8 x 10" 8 . The resultant Kerr 
angle (j2~4")l is 6k = 230 nanorad. This estimate is 3.6 
times greater than the experimentally observed value of 
65 nanorad [l[. The experimental Kerr angle may be re- 
duced relative to the theoretical estimate for variety of 
reasons. For example, the effective value of Ao at high 
energies may be lower than at the Fermi level. We con- 
clude that the theoretical formula (f2~4"|) reasonably agrees 
with the experiment by the order of magnitude. 

A different theoretical formula with Ok oc Ao was pro- 
posed phenomenologically in Ref. [l| motivated by the 
experimental temperature dependence of 0k(T). On the 
other hand, our formula (f2~4")l gives Ok oc Aq. The error 
bars in the experiment pjare quite big, so deciding be- 
tween the linear or quadratic dependences of Ok on Ao 
may require more precise measurements. The appearance 
of Aq in Eqs. (|2"Tj) and (|24|) is quite natural, originating 
from the product A x A y , which changes sign when the 
chirality of the order parameter changes from p x + ip y 
to p x — ip y , as observed experimentally 1]. This product 
can be also written as the vector fx** pointing along z 
14j . which is consistent with Eq. (fT7| . It would be very 
interesting to verify experimentally the Q~ 3 frequency 
dependence of the Kerr angle predicted by Eq. (j2"4")) . 

Experiments indicate that the superconducting gap 
may have the so-called horizontal lines of nodes in 
S^RuCU (see Ref. jll| and references therein). In this 



case, Ao should be considered a function of the elec- 
tron momentum p z perpendicular to the layers: Ao — > 
A cos(p z d) for tri plet pairing [ll[ or A — * Aosin(p z d) 
for singlet pairing [lfj. In both cases, averaging A^(p z ) 
over p z generates an additional factor 1/2 in Eq. (|2"4|) 
and no changes in Eq. (fT9|) . Thus, although experi- 
ment [l| directly proves the chiral character of the su- 
perconducting pairing in Sr2Ru04, it does not discrimi- 
nate between triplet pairing and the chiral singlet pairing 



(j>x + iPy) sva.{p z d) proposed in Ref. [loj ]. 

In conclusion, we derived the Chern-Simons-like term 
in the effective action of the 2D chiral p x + ip v super- 
conductor, generalizing previous results @, [IJ to finite 
frequency and temperature. The resultant dc Hall con- 
ductivity has the half-quantum value a xy — e 2 /2h at 
T = [6j, but is reduced at a finite temperature by the 
factor (|2U)l . We derived Eq. (f2"Tj) for the imaginary part 
of the optical Hall conductivity and Eq. (|24|) for the po- 
lar Kerr angle, which agrees by the order of magnitude 
with the recent experimental measurement in Sr 2 Ru04 
by Xia et al. [l|. Eq. predicts that the Kerr angle is 
proportional to the square of the superconducting energy 
gap and is inversely proportional to the cube of frequency, 
which can be verified experimentally. The derivation may 
be also relevant for the finite-temperature Chern-Simons 
theories in high-energy physics (see references in 2(J 2lj|). 
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